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Abstract - We study the joint statistics of conductance G and shot noise P in chaotic cavities 
supporting a large number N of open electronic channels in the two attached leads. We determine 
the full phase diagram in the ( G , P) plane, employing a Coulomb gas technique on the joint density 
of transmission eigenvalues, as dictated by Random Matrix Theory. We find that in the region of 
typical fluctuations, conductance and shot noise are uncorrelated and jointly Gaussian, and away 
from it they fluctuate according to a different joint rate function in each phase of the (G, P) plane. 
Different functional forms of the rate function in different regions emerge as a direct consequence 
of third order phase transitions in the associated Coulomb gas problem. 


cd 


Introduction. — We consider the problem of quan- 
q turn electronic transport in mesoscopic devices. The typ- 
Qical setting is a cavity etched in semiconductors and con- 
L ““ J nected to the external world by two attached leads. The 
cavity is brought out of equilibrium by an applied exter- 
nal voltage. Assuming that the average electron dwell time 
(spent inside the cavity) is well in excess of the Ehrenfest 
\| time (when quantum effects kick in) at low temperatures 
^ and applied voltage, statistical properties of the electronic 
, V. transport through a cavity exhibiting chaotic classical dy- 
• namics display remarkably universal features. 

Random matrix theory (RMT) has provided tools and 
invaluable insights to successfully describe this universal- 
t-H ity 0- In the RMT approach, the scattering process inside 
l ' the cavity is governed by a scattering matrix S drawn at 
• i—i random from the unitary group [2], 3 . 

Vn The two incoming leads may in general support N\ and 
V 2 (with Ni < IV 2 ) open electronic channels, i.e. differ¬ 


ed 


The S matrix, having a natural block form S = 


ent wave numbers of the incoming electronic plane waves. 

J' r '. 

in terms of reflection (r, r') and transmission [t,t') matri¬ 
ces, connects the incoming and outgoing electronic wave- 
functions. Conservation of electronic current implies that 
S is unitary. Landauer-Biittiker theory [4jj6] expresses 
most physical observables in terms of the eigenvalues 


A = (Ai,..., AjVi) of the Hermitian matrix tt'. Unitar- 
ity of S implies that 0 < A,; < 1, V*. Hereafter, we will 
consider symmetric cavities N = N\ = N- 2 . 

Important observables are the (specific, per channel) 
conductance and shot noise 

P=^5>(l-A,), (1) 

i i 

measured in units of the conductance quantum Go = 
2e 2 //i and of Po = 2e|AF|Go (with AV the applied volt¬ 
age), respectively. They are expressed as linear statis¬ 
tics on the transmission eigenvalues. Since the eigenval¬ 
ues are A^ = 0(1), both G and P are intensive quantities 
of order 0(1). More precisely, 0 < G < 1 and, since 
Ei A ? < (£* Ai) 2 < N )Tk A 2 , we have the useful geomet¬ 
ric inequalities 


0 < P < G(1 — G) < 1/4. (2) 

The upper bound is attained at G = 1/2, while the lower 
bound at G = 0 or G = 1. 

Assuming now that S' is a random unitary matrix then 
implies that the AjS become (correlated) random variables. 
What is their distribution? When the leads attached to 
the cavity are ideal (no tunnel barriers), the appropri¬ 
ate choice is to assume S uniformly distributed in one of 
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Dyson’s circular ensembles of random matrices, labeled 
by a parameter j3: S is unitary and symmetric for p = 1 
(for systems that are invariant under time-reversal), just 
unitary for p = 2 (broken time-reversal invariance) and 
unitary self-dual for p = 4 (in case of anti-unitary time- 
reversal invariance). In this ideal case the joint probability 
density function (jpdf) of transmission eigenvalues is given 
[2,7,8] by the Jacobi ensemble of RMT, namely 

/(A) = ^ni^-A/lK /2 -\ A €[0,1]*, (3) 

1<J 2=1 


where Zn is a normalization constant enforcing 
Jj 0 dA/(A) = 1. The random nature of A promotes 
conductance and shot noise to random variables them¬ 
selves, whose statistics is of paramount interest. The aver¬ 
age and variance of conductance and shot noise were con¬ 
sidered, using perturbation theory in 1 /N, long ago mu 
In particular, as N —> oo the variance does not scale with 
IV, as one would naively expects, but instead attains a 
constant value (oc 1/P) depending only on the symmetry 
class, a phenomenon that has been dubbed universal con¬ 
ductance fluctuations. This value can be predicted from 
classical variance formulas 10 . More recently, the classi¬ 


cal theory of Selberg integral was employed and extended 
ll|jl4] to address the calculation of transport statistics 
non-perturbatively (i.e. for a fixed and finite number of 
channels). The full distribution of G and P is strongly 
non-Gaussian for few open channels, with power-law tails 
at the edge of their supports and non-analytic points in 
the interior 14-16 . For finite N and p = 2, the Laplace 


transform of the full distribution was studied in 17 spot¬ 
ting a connection with integrable systems and Painleve 
transcendents. The full distribution (including large devi¬ 
ation tails) for large N was studied in 18] using a Coulomb 
gas technique, where an error occurring in the asymptotic 
analysis in [l7j was corrected. The statistics of other ob¬ 
servables was studied in 13[l9-22 . The integrable theory 


of quantum transport in the ideal case, pioneered in 17 


for P = 2 , has been recently completed including the other 


symmetry classes 23 


In this paper, we are concerned with the joint statistics 
of conductance and shot noise 0 for a large number of 
open channels N and in both regimes of small (typical) 
and large (atypical) fluctuations. Of particular interest 
are the jpdf of G and P 


V(g,p) = (5(g-G) S(p-P)) (4) 


(where (•) stands for the average with respect to the 
jpdf and its Laplace transform 

V(s,w) = (e-^ sG+wP ^. (5) 

As discussed above, the arsenal of sophisticated tech¬ 
niques employed to tackle the calculation of joint cumu- 
lants for finite N, and their leading asymptotic behavior 


for N —> oo, is truly impressive and virtually leaves no 
room for improvement in terms of mathematical rigor. In 
physical terms, however, it is desirable to get a more in¬ 
tuitive understanding, based on simpler and more imme¬ 
diately decipherable formulas. Here we use the physically 
transparent Coulomb gas method to get (at least to our 
eyes) a neater picture of the mutual relations between con¬ 
ductance and shot noise in such systems. The obtained 
formulas for the large deviation functions allow to recover 
effortlessly the leading terms of the joint cumulants. Ad¬ 
ditionally, the typical fluctuations of functions of both G 
and P can be easily investigated. As an example we will 
provide results on the Fano factor 

f = 5 < 6 > 


which is essentially the ratio of the actual shot noise and 
the Poisson noise that would be measured if the system 
produced noise due to single independent electrons (for 
more details see the review 24 ). Moreover, the full phase 
diagram in the (G, P) plane is obtained exactly and linked 
to the equilibrium electrostatic properties of the associ¬ 
ated Coulomb gas of charged particles. We begin by first 
summarizing our results. 


Summary of results. — In this work we show, using 
a Coulomb gas technique that, for large N and any p > 0, 
the Laplace transform ([5| behaves as 

V(s,w) ( 7 ) 


where J(s, w), independent of p, is the generating function 
(GF) of the joint cumulants of G and P. Henceforth, the 
symbol ss stands for equivalence on a logarithmic scale. 
From this result, large deviation theory predicts that the 
jpdf of (G, P) in Q behaves asymptotically as 

P(g,p)ne-% N2 *^\ ( 8 ) 


where T( g,p) = inf ajW [J(s, w) — sx — wp\ is the so-called 
joint rate function. We find that J(s, w) (and correspond¬ 
ingly *F(g,p)) takes five different functional forms in differ¬ 
ent regions of the (s, w) (respectively, ( g,p )) plane. This is 
a direct consequence of phase transitions in the associated 
Coulomb gas problem (see next Section). Across the lines 
of phase separation, the third derivatives of the GF (the 
free energy of the associated Coulomb gas) are discontin¬ 
uous, as it is typical in this type of problems (see [25] for 
a recent review). The rate function ’F( g,p ) has a global 
minimum (zero) at (g,p) = ( 1 / 2 , 1 / 8 ), corresponding to 
the average value of conductance and shot noise for large 
N, (G) = 1/2 and (P) = 1/8. Expanding the rate func¬ 
tion around this minimum, we find that the typical joint 
fluctuations of conductance and shot noise are Gaussian, 
with a diagonal covariance matrix [see Eq. (25)] implying 
the absence of cross-correlations to leading order. This 


is in agreement with earlier findings 12 14 17 19,23 26 
However, atypical fluctuations far from the average are 


P-2 

























Joint statistics of quantum transport in chaotic cavities 


not described by the Gaussian law, but rather by a dif¬ 
ferent rate function. For w = 0 (s = 0) our GF reduces 
to the GF of conductance (shot noise) alone, computed 
in 18 . In the next Section, we set up the Coulomb gas 


Originally due to Dyson 28 , this Coulomb gas tech¬ 


nique with additional constraints has been developed and 
used in many different problems 18p5||29 - 36 which have 


calculation. 

The Coulomb gas. — The suitable tool to perform 
a large dimensional analysis is the empirical density of 
transmission eigenvalues 


p(A) = iV- 1 ^5(A-A i ) 


(9) 


apparently very little to do with each other. Recently, an 
improvement of the method has been used to compute 
the joint statistics of two linear statistics in the Wishart- 
Laguerre ensemble of random matrices 37 38 . We will 
adopt the same strategy below. 

First, it is worth noticing a symmetry of this problem. 
The external potential is invariant (up to an immaterial 
constant) under the exchange (s, A) —> (—s, 1 —A), namely, 


a normalized random measure on the interval [0,1]. The 
first step to compute the multiple integral ([5]) is to realize 
that it can be written as T{s,w) = w)/Zn(0, 0), 

where 

( 10 ) 


T4xt(l - A ;-s,w) = I4xt(A ;s,u>) - s. 


(15) 


Z N (s,w) = J [Dp] 

and the energy density functional 


e -f N 2 E[p-,s,w] 


E[p\ s,w] = — Jj log |A — A'lp^/^A^dAdA' 
+ J 14 xt(A; s, w)p(X)d\ + 0(l/N) , 


Since the logarithmic repulsion is also invariant under the 
exchange A —> 1 — A, the energy functional E[p;s,w\ in¬ 
herits this symmetry. Therefore, the equilibrium spectral 
densities at s and — s are simply related by p*( A; —s, w ) = 
p*( 1 — A ;s,w). Thus the phase diagram is invariant un¬ 
der the inversion s —> — s and it is sufficient to study the 
problem for s > 0. In the next Section, we compute the 


explicit solutions of the singular integral equation (141 for 


( 11 ) 


any given value of ( s,w ) and we find that five different 
functional forms are possible for the saddle-point density 
p* in various regions of the (s, w) plane. 

with 14 xt (A;s,w) = (s + in) A — w\~, is^defined onmormal- Phase diagram. — We have identified five domains 

flj, 1 < j < 5 in the (s, u>)-plane that correspond to five 
different phases of the Coulomb gas. The functional forms 
of p* and their supports in the phase ( s,w ) € flj , with 


ized spectral densities / p(A)dA = 1. Equation (10) is the 
Gibbs-Boltzmann partition function of a system of charged 
particles confined on the interval [0,1] of the real line, in¬ 
teracting via a (logarithmic) 2D Coulomb potential in an 
external single-particle potential 14xt ■ The long-range, all- 
to-all nature of the interaction between the charges has 
the consequence that both energy contributions (interac¬ 
tion and external potential) are of order 0(N 2 ). This is 
the origin of the 0{N 2 ) speed of the large deviation func¬ 
tions Q and in sharp contrast to the standard O(N) 
speed for classical large deviation theory of independent 
random variables 27 . Using p, any linear statistics for 


j = 1,..., 5, read 


Pi(A) = 
P2W = 


91(A) b- A 


2i V A 

92(A) 1 

2 tt ^/A(l — A) ’ 
93(A) 


A €[0,6], (16) 

A €[0,1], (17) 


P 3 ( A) = ^\AA-«)(&-A), A € [a, 6 ], (18) 


large N can be easily written down. In particular, 

G\p] = J Ap(A)dA, P[p\ = J A(1 - A) p(A)dA . (12) 


pi(A) = MMj (a * )(A b \ A e [0,a] U [6,1], (19) 


PU A) = 


In the large N limit, the partition function (101 is dom¬ 
inated by the saddle-point density 


2tt y A (1 - A) 
95(A) jX - a 

2 t r V 1 - A’ 


A € [a, 1], (20) 


Z N (s,w) ~ e ~ 2 N2E lp*'’ s ’ w ]. 


(13) 


with 0 < a < b < 1. Here 9 j(A) are polynomials in A, 
whose expressions, together with those of a and b , depend 


Here p*( A) (depending parametrically on s and w ) is the 
solution of the saddle point equation j^E[p\s,w\ = 0, 
namely 

± Vext (\;s,w) = 2j : (14) 

for all A where the charge density exists p*{ A) > 0, where 
f denotes Cauchy’s principal value. In terms of our elec¬ 
trostatic model, Eq. (14) is the continuous version of the 


force balance condition for the charge cloud p*( A) to be in 
equilibrium. 


on the values (s,w) and are given explicitly in (31)-(32) 
of the Appendix. The five different domains and the cor¬ 
responding densities are shown in Fig. [T] 

The physical picture is quite intuitive in terms of the 
interplay between the logarithmic interaction between the 
charges and the external potential U e xt- For moderate 
values of s and w (i.e. (s,w) € U 2 ) the external potential 
is too weak compared to the logarithmic repulsion between 
the charges, so that the gas invades the whole interval 
[0,1]. For sufficiently large positive values of s (Hi), the 
potential is strong enough to attract the charges towards 
A = 0 and confine the gas in [0, b] C [0,1]. In a similar way, 
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Fig. 1: Top Left: Phase diagram in the Laplace ( s , w ) plane. The five domains fij, i = 1, ..., 5 correspond to different spectral 
profiles p*(A) of the transmission eigenvalues (|16|)-(|20|). The critical lines that separate the phases are reported in Appendix 


in (281-( 301 . Notice that the phases fb (suppp* = [0,6]) and fls (suppp* = [a, 1]) are separated by the phases flo, S ?3 and 1 ?4. 
The diagram is symmetric with respect to the line s = 0. The lines w = 0 (black dotdashed) and s = 0 (gray dashed) have been 
studied in £[8]. We recover their critical points s cr = ±4 and ui cr = ±8. Top Right: The corresponding phase diagram in the 
conductance/shot noise (g,p)-plane. For completeness, the critical points and the lines corresponding to s = 0, w = 0 are also 
reported in real space. Bottom: Different spectral profiles p*( A) of the Coulomb gas. The analytical curves are superimposed 
to points from a Montecarlo simulation of a Coulomb gas of N = 30 particles in the external potential 14 x t(A; s, w). 


large negative values of s (SI 5 ) provide a repulsive potential 
that pushes the gas away from the origin towards the hard 
wall located at A = 1. When w < 0 and |s| < — w, I4xt 
has a minimum that attracts the charges; if this well is 
sufficiently deep (i.e. (s, w) £ V2 3 ) the Coulomb gas is 
trapped in the minimum of the potential and does not 
feel the hard walls. On the other hand, when w > 0 and 
|s| < — w , I 4 xt has a maximum that repels the charges; in 
1^4 the peak is sufficiently high to cause the Coulomb gas 
to split into two separated components that tend to stay 
far apart. 


The single-support solutions have been found employ¬ 
ing a theorem due to Tricomi 39 for the singular equa¬ 
tion ( fl4| . The theorem provides the general form of 
p*( A) = p*(A;C, a, 6 ) with a single support. It depends 
on three arbitrary constants that are determined by im¬ 
posing the normalization condition C = f p*( A)dA = 1 
and the behavior of p* at the two edges a and 6 of the 
support. This procedure provides a density supported on 
a single interval and parametrized by s and w. 

Multiple-support solutions of the saddle-point equations 
are more complicate and involve in general more arbi¬ 
trary constants. For example for a solution p*(A) = 
p(A; C, a, b, n) with double support in [0, a] U [ b , 1] (the 
only multiple support scenario that we have in our prob¬ 


lem), we have to impose the overall normalization C = 1, 
the regularity condition at the edges p*(a) = p*(b) = 0 and 
an extra constraint, i.e. the so-called filling fraction of one 
of the two component of the density, ff p*(A)dA = n 
How to fix this filling fraction? At equilibrium, the 


ing fraction k is such that a small variation 5n does not 
change the energy at first order. This condition means 
that, at equilibrium, the work required to move a single 
charge from a to b is null, and can be easily pro ved to reac 0 
ip*(A)dA = 0 where suppp* = [0, a] U [b, 1] 


40 


41 


These kinds of constraints are often cumbersome to eval¬ 
uate. However, if one is only interested in finding the 
critical line of phase transition between the single-support 
phases and the double-support phase in the (s, ru)-plane, 
a practical strategy is the following. The condition of null 
work for moving a charge from a to b should be satisfied 
for any double-support solution. In the limit this will be 
valid also at the “birth” of the double-support phase when 
6=1. Then, the requirement J p*(A)dA = 0 with p*{ A) 
supported on [ 0 , 6 *] provides a threshold 6 * and therefore 
an equation connecting s and w at the “birth-of-second- 
cut”, i.e. precisely the line of phase transition in the (s, w)- 


1 Note that the integral runs outside the support of p*, therefore 
it needs to be interpreted as the integral of the analytic continuation 
of p*(A) to values of A not in the support. 
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plane. 

For given values of ( s,w ), the solution p* provides the 
typical configuration of eigenvalues yielding prescribed 
values for G and P from (12). Both are functions of s 
and w: 


joint cumulant GF J(s,w) follows by a careful integration 
of the differential form ( 22 ) and the continuity require¬ 
ment through the critical lines. By starting from the cen¬ 
ter ( 0 , 0 ) of region f l 2 , where J( 0 , 0 ) = 0 , and integrating, 
one obtains 


g{s,w) = G[p *], p(s,w) = P[p*]. 


( 21 ) 


In particular, for (s, w) = (0, 0) (uncostrained problem) 
we obtain from 0 and the arcsine law p as (A) = 
l/^y/ACl — A)) providing the average values ( G ) = 1/2 
and (P) = 1/8 from ( fl2| ). In general, © in Go describes 
the typical fluctuations around the average values ( G ) and 
(P). Values of (s,w) £ Hi correspond to a configuration 


of the transmission eigenvalues (16) yielding smaller val¬ 


ues of the conductance g (G) while large deviations 
g ( G) correspond to (s,w) £ Sis and are described 

by (p0|. Similarly, large values for P are given by a typ¬ 


ical configuration of eigenvalues as in (18) in phase H 3 , 


while smaller values for P in phase fl 4 are ascribed to the 


double support solution (19). As the point (s, w) moves in 
the Laplace plane, the configuration of eigenvalues p*{ A) 


changes according to (16)-(20). The transition across the 


critical line 7 ^ (the boundary of two phases f \ and Clj) 
corresponds to a change in shape of the Coulomb gas den¬ 
sity, which signals a corresponding phase transition in the 
quantities G and P. In the next Section, we use the ex¬ 
plicit functional forms of the density in various regions of 
the (s,w) plane to show how the GF and the rate functions 
can be computed. 

Joint large deviation functions. First, one in¬ 
serts the equilibrium density p* (whose explicit expres¬ 


sions are given in Eqs. (16)-(201) into (12) and computes 
the corresponding single integrals, yielding the order 0(1) 
quantities g(s,w) and p(s,w). These are in turn related 
to the cumulant GF J(s,w) in ([ 7 ]) and the rate function 
d>(g,p) in (|8]) via the differential relations 

dJ(s, w) = g(s, w)ds + p(s, w)dw , (22) 

-d>F(g,p) = s{g 1 p)dg + w{g,p)dp , (23) 

complemented with the condition J(0,0) = 0 and 

^( 3 ( 0 , 0),p(0,0)) = JF(1/2,1/8) = 0. Here s(g,p) and 
w(g,p) are the solution of s = g(s,w) and w = p(s,w). 
The expressions (22) and (231 are known as Maxwell re¬ 


lations among thermodynamic potentials, in our case the 
Helmholtz free energy and the enthalpy. The application 
of standard thermodynamics arguments to the Coulomb 
gas can be rigorously justified by means of large deviation 


principles 42 . In fact, relations (22) and (23) are equiva¬ 


lent statements: one is in Laplace space - the other one is 


in real space. By plugging (16)-(20) into (21), one easily 


obtains the explicit expressions of g(s,w ) and p(s,w) in 
each region of the phase space. The lines w = 0 and s = 0 
(yielding the GF of conductance or shot noise alone) were 
studied in [18] and we easily recover the known results by 
taking the limit w —> 0 (or s —» 0) in our expression. The 


+ (24) 

Therefore, in the region of typical fluctuations, G and 
P are (to leading order) uncorrelated Gaussian variables 
with average values (G) = 1/2 and ( P ) = 1/8, VarG = 
l/(8/37V 2 ), VarP = 1/(64/3iV 2 ), and Cov(G,P) = 0: 

P(G, P) « g-4£JV 2 (G-|) 2 -32/3iV 2 (F-|) 2 _ (25) 


Since the higher derivatives of J(s, w ) vanish at (s,w) = 
(0, 0) we also conclude that the higher order cumulants 
are P) = G(l/7V £+m - 1 ) for l + m > 2. 

For generic s and w, the integration of dJ provides the 
cumulant GF: 

r(s,w) 

J(s,w) = J(s 0 ,w 0 ) + / dJ(s / , w'), (26) 

J (sqj'U’o) 


for (s,w) £ fij. Here, one has to choose a suitable ini¬ 
tial point on the boundary with a phase f tj already com¬ 
puted, namely (s o,wq) 


£ 


D f tj. The expression (261 


can be evaluated by integrating in the complex plane and 
applying the residue theorem. Across the lines of phase 
separation, one finds that the GF is not analytic, since its 
third derivatives are discontinuous. 

As an additional bonus, from the joint limiting behav¬ 
ior of G and P we can deduce the distribution of functions 
of both G and P like the Fano factor §. We stress the 
fact that F is not a linear statistics of the transmission 


eigenvalues. From (25) we conclude that the typical fluc¬ 


tuations of F are asymptotically Gaussian 

P(P)«e-^ Jv2 ( F -^) 2 , 


(27) 


with average value ( F) = 1/4 and and variance VarP = 
3/(32/3V 2 ). To the best of our knowledge the asymptotic 
distribution ( [27] ) is a new result. 

Conclusions. — In summary, our analysis gives an 
overall picture of the joint statistics of conductance G and 
shot noise P for an ideal chaotic cavity supporting a large 
number N of electronic channels in the two attached leads. 
We employed a Coulomb gas technique to establish the 
large deviation formulas ([T]) and ([8]) , governing the behav¬ 
ior of the joint cumulant GF and the joint rate function of 
G and P. We were able to obtain the full phase diagram 
in the (s, w) and (G, P) planes and we found that both 
J(s,w) and d<(g,p) acquire five different functional forms 
in different regions of their domain. These different ex¬ 
pressions are a direct consequence of phase transitions in 
the associated Coulomb gas problem. Across the lines of 
phase separation, the third derivatives are discontinuous, 
implying a third order phase transition. 
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Appendix. — Let us denote by 7 *j the critical curve 
that separates the phases fit and f lj in the (s, w) plane. 
We report their explicit expressions 

713 U 735 : |s| = —w — \/— 8 w for w < —8 , 

712 U 725 : |s| = w/2 + 4 for — 8 < w < 8/3 , 

724 : s 2 + 2w 2 — 16w = 0 for 8/3 < w < 8 , (28) 

while the curves 714 U 745 are the solutions of the equation 

2wcj) 1 ( 6 !(|s|, tu)) = [2(|s| + w) + \J (|s| + w) 2 - 24u>] 

X </>0 (h{\s\, w)) /3 , w> 4, (29) 

(f>k{b) = i f &xx k \f (b — x)/x , 0 < b < 1 . (30) 

Jb 

The polynomials < 7 ,; (A) and the edges a, . b, of the sup¬ 
ports suppp* are given by 

91 (A) = [- 6 tuA + 2 (s + w) + y/(s + w ) 2 - 24uj]/3 , 
92 (A) = 2wA~ — (2 w -f- s)A T wj '4 -t- s/2 T 2 , 

93 (A) = - 2 w , 

( 74 (A) = 2 (t — wx) , 

95 (A) = 7 i (1 — A; ~s, w) , (31) 

and 

cii = 0, bi = [s + w — \/(s + w) 2 — 2Aw]/(3w), 

02 = 0 , &2 = 1 , 

a 3 = (s + w + V—8w)/(2w), 63 = 1 - a 3 (—s,w), 

04 = [s + (2 — 64 ) 11 ; — 2 t]/w, 

64 : / 2 (t — wx)\J (04 — x)(x — 64 )/\/x (1 — x)dx = 0 , 

J CI4 

t = (s + 2w — b^w) /3 

— \Js 2 + 2(264 — l)sw — w (24 — w — 86411 ; + 8b1w)/6, 
a 5 = l- bi(-8,w), 65 = 1 , (32) 

respectively. 
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